Probing the Binary Black Hole Merger Regime with Scalar Perturbations 
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We present results obtained by scattering a scalar field off the curved background of a coalescing 
binary black hole system. A massless scalar field is evolved on a set of fixed backgrounds, each 
provided by a spatial hypersurface generated numerically during a binary black hole merger. We 
show that the scalar field scattered from the merger region exhibits quasinormal ringing once a 
common apparent horizon surrounds the two black holes. This occurs earlier than the onset of 
the perturbative regime as measured by the start of the quasinormal ringing in the gravitational 
waveforms. We also use the scalar quasinormal frequencies to associate a mass and a spin with 
each hypersurface, and observe the compatibility of this measure with the horizon mass and spin 
computed from the dynamical horizon framework. 
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I. INTRODUCTION 

Over the past thirty years, numerical evolutions of bi- 
nary black hole (BBH) systems have provided deeper and 
deeper insight into the properties of the vacuum two- 
body problem ofgeneral relativity. Today a number of 
codes [llSSHjHIESIl] are able to accurately evolve 
black hole binaries; in particular, the merger regime it- 
self is now available for studies. It was shown in refer- 
ence [13] that the merged system seems to enter the 
perturbative stage while it is still evolving and radiat- 
ing energy and angular momentum. This is in agree- 
ment with earlier results from the Close Limit approxi- 
mation [ll| and from the Lazarus project p^ . 

In order to study the transition between the nonlinear 
merger stage and the final ringing black hole, we scatter 
a test scalar field off the geometry generated in the in- 
spiral and merger of the binary system. Test fields and 
perturbation theory have traditionally been used to study 
excitations on a single stationary black hole spacetime, 
both for the Schwarzschild [H, [3, [3 and for the Kerr 
solutions p^ . Numerical methods were used in several 
scenarios (see for instance [13] and ^18]); more recently, 
these results have been extended to more general black 
hole solutions, including non- vacuum spacetimes [l9l| and 
models with non-flat asymptotic topology [13, l2ll|. 

An interesting result of quasinormal mode (QNM) evo- 
lution on non- stationary, spherically symmetric black 
hole spacetimes is that the mass extracted from the com- 
plex quasinormal frequency tracks the time-dependent 
black hole (BH) mass quite closely (modulo a delay ef- 
fect), equaling the constant quasinormal frequency one 
would obtain on a static spacetime with the mass param- 
eter equal to the instantaneous mass at each time [l9| . 
This adiabatic picture has an intriguing parallel in black 
hole thermodynamics: the application of the first law to a 
spacetime containing an axisymmetric dynamical horizon 
shows that the mass increase between two cross sections 
with horizon area Ai and A2 and angular momentum 



Ji and J2 is given exactly by Af(A2, J2) — M{Ai,Ji), 
where M{A, J) represents the functional dependence of 
the mass M on A and J for a Kerr solution [22|. In 
other words, as far as the first law of thermodynamics 
is concerned, this portion of the spacetime can formally 
be pictured as a stack of spatial hypersurfaces, with each 
slice belonging to a member of the Kerr family with a 
different mass and spin parameter. 

In this study we aim to address some questions con- 
cerning the transition between a system of two inspiraling 
black holes and a single perturbed black hole. We inves- 
tigate how early during the BBH evolution the scalar 
field probe exhibits the quasinormal ringing phase or- 
dinarily expected from a spacetime containing a single 
black hole, and what features this phase displays when 
compared to gravitational wave quasinormal ringing of 
the single stationary black hole that is left behind by the 
coalescence process. Based on these considerations we 
present a novel strategy that consists of (i) simulating 
a BBH merger and extracting horizon information from 
the merged black hole, (ii) taking snapshots of the spatial 
geometry obtained in the BBH merger at a sequence of 
coordinate times, (iii) evolving a massless scalar field on 
the corresponding spatial hypersurfaces, and (iv) observ- 
ing the scalar field undergo a quasinormal ringing phase 
and extracting its complex frequency, which in turn can 
be converted into a mass and a spin parameter using the 
functional dependence of the quasinormal frequencies on 
these two parameters prescribed by perturbation theory. 
We then compare these estimates to the mass and spin 
extracted directly from the horizon. 

This method has the benefit of effectively separating 
the evolution time of the binary system from the relax- 
ation time of the scalar perturbations, allowing for a mea- 
sure of the "instantaneous" properties of the foliation. 
Notice that this procedure is effectively equivalent to the 
evolution of a scalar field on the fictitious spacetime one 
would obtain by selecting each spatial surface and freez- 
ing it in time. This of course is not meant to constitute a 
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physical solution to the coupled scalar-gravitational sys- 
tem: as emphasized above the scalar field is used here 
exclusively as a probing agent for individual hypersur- 
faces. 

We will present a brief overview of scalar perturbations 
in section ini along with the method we use to evolve the 
scalar wave equation on a numerically-generated curved 
background, extract the complex frequencies and asso- 
ciate a mass and spin parameter to them. Section IIIII 
contains the results of our scalar field evolutions, along 
with the corresponding waveforms and analysis of errors. 
Finally, in our concluding section Hvl wc will discuss the 
results. 



II. SCALAR FIELD EVOLUTION 

The behavior of a massless scalar field on the back- 
ground of a single Kerr black hole has been the subject 
of extensive analytical treatments. A good introduction 
is given in reference [IS] and a thorough treatment of the 
subject can be found in NoUert's |2^] and Kokkotas and 
Schmidt's review works 12^. 



A. Quasinormal ringing 

The behavior of a massless scalar field $ on the back- 
ground provided by a Kerr black hole of mass M and Kerr 
spin parameter a = J/M is governed by the source-free 
scalar wave equation: 



□$ := gabVV''<S> = 0. 



(1) 



Separating out the angular variables, one reaches a wave 
equation with a potential [Id^ . Upon further study, one 
finds that the field exhibits many of the well-known fea- 
tures of waves impinging on a potential; for instance, it 
exhibits quasinormal ringing, i.e. a phase of decaying os- 
cillations characterized by a certain complex frequency 
spectrum. 

The complex nature of the frequency spectrum leads to 
several undesirable consequences, some of which have di- 
rect ramifications on our ability to detect and character- 
ize quasinormal ringing in a numerical simulation: first 
and foremost, the spectrum of quasinormal modes does 
not constitute a complete set of eigenfunctions [H, [H, 
[26j : early on, the behavior of the scalar field depends on 
the process that generated the perturbation, and the field 
is said to undergo the prompt response regime, where $ 
still displays memory of the perturbing mechanism [131 • 
Additionally, at late times the backscattering off the cur- 
vature will overtake the by-now feeble damped oscilla- 
tions and dominate the waveforms, generating the so- 
called tails 27[. A Green's-function analysis of the ori- 



question then arises about how to identify the onset time 
of quasinormal ringing, since an accurate determination 
of this quantity is of the utmost importance while fitting 
for the quasinormal frequencies. This issue is usually re- 
ferred to as the time shift problem; for a discussion of its 
influence on the bias of the fitted parameters, see p^ . 
Also, concepts such as the fraction of a given waveform 
that can be identified as each quasinormal mode and the 
gravitational energy contained in each quasinormal mode 
cannot be defined in rigorous terms |30|; the absence of 
such auxiliary notions also poses serious complications to 
identifying the onset of quasinormal ringing in a numer- 
ical waveform, and quantifying its quasinormal content. 
Berti et al. [1^ review this issue from the numerical point 
of view. 

Notice that, due to their dependence on the mass M 
and spin parameter j = a/M oi the background Kerr so- 
lution, the quasinormal ringing portion of the scattered 
waves carries information about M and j away from 
the scattering center (a property which, for instance, 
makes gravitational quasinormal ringing a powerful tool 
for gravitational wave astronomy 31]). 

In this work, we use scalar perturbations to probe the 
parameters describing the final black hole that results 
from the binary coalescence. In order to provide as ac- 
curate a depiction of the black hole as possible, we will 
first have to devise a strategy to tackle the ambiguities 
associated with quasinormal modes and find a reliable 
strategy to detect and measure the quasinormal compo- 
nents in a waveform. In section III El we will present our 
experimental method. 



B. Scalar field evolution on a numerically 
generated background 

We evolve a scalar field on a numerically generated 
background obtained from a binary black hole evolution. 
We take snapshots of the spatial geometry and then use 
each of these as a time-independent background on which 
the scalar field propagates. 

During the course of the evolution, we store the values 
of the grid functions corresponding to the gravitational 
variables at a sequence of coordinate times {tAr}, which 
we denote with a bar: 



K,j{x'') = Kij{tN,x'') 
a{x'') = a{tN,x^) 
P\x'')=(3HtN,x'') 



(2) 
(3) 
(4) 
(5) 



gin of these three phases, along with a convergence study 
of the quasinormal mode exp ansion is presented in (28| 
for Schwarzschild and in [2^ for Kerr black holes. The 



and evolve a scalar field on the background provided by 
the corresponding frozen geometry. In this sense, $ is 
nothing more than a probing agent for each of the spa- 
tial slices on which we perform the experiment, and is 
meant solely to provide information, through its scatter- 
ing, about the curvature of each single slice rather than 
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to explore a scalar-gravitational coupling of an actual 
physical system. 



C. Background spacetime 

The evolution of the gravitational variables Tjw and Kj^ 
is performed according to the BSSN system [sl, [H, Hj] , 
as implemented in the PSU code through the Kranc gen- 
erator [sHi- The fourth order finite differencing code is 
based on the Cactus infrastructure [S^ and is endowed 
with Carpet's mesh refinement capabilities [37|. The 
gauge variables a and are initialized and evolved ac- 
cording to the moving puncture prescription 0, Q , while 
the initial data are generated following [sst - 

For the binary coalescence runs, we choose an initial 
configuration usually referred to as Rl [s^: this setup 
represents two non-spinning black holes with irreducible 
masses equal to 0.505A1 (where A4 denotes a mass scale 
related to the ADM mass via Eabu/M = 0.9957). 
The two black holes start at a coordinate separation of 
6.514A^, with linear momenta perpendicular to the sep- 
aration vector and equal to 0.133M. Further details re- 
garding our code and our earlier study of this system can 
be found in We used a reference spatial resolution 
of A^/51.2, and evolved some simulations at two extra 
resolutions (A1/44.8 and Ai/57.6) to obtain an estimate 
for the truncation errors (see section UlI CI below) . 

Figure [T] plots the gravitational wave from the simu- 
lated system: the two black holes inspiral around each 
other for about two orbits before plunging and forming a 
single final perturbed black hole whose parameters can be 
inferred, due to energy and angular momentum conserva- 
tion, from the emitted radiation in terms of the Newman- 
Penrose scalar \['4, as discussed in [41|]. For the M/51.2 
resolution run, starting with an initial ADM mass and 
angular momentum equal to Eabm/M = 0.9957 and 
Jabm/M'^ = 0.862, the system radiates E^i^^/M = 
0.032 ± 0.005 and Jrad/X^ = 0.22 ± 0.01, leaving be- 
hind a final black hole with Ef/M = 0.963 ± 0.005 and 
Jf/M.^ — 0.64±0.01, in agreement with the estimates of 
mass and spin provided by the analysis of the fundamen- 
tal quasinormal tone of "^f (fi'QNivi/A^ = 0.957 ± 0.005 
and Jqnm/A^^ = 0.62 ± 0.01) and with previous results 
obtained with other codes [391 . 



D. Scalar evolution in 3+1 form 

The scalar field code was developed using Kranc [sH] 
to generate the initial data and evolution modules and 
the interface with the BSSN code. A 3-1-1 decomposition 
of the scalar wave evolution equation ([1]) in terms of 7^ , 
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FIG. 1: The real part of the ^ = 2, m = 2 mode of the 
Newman-Penrose scalar 'I'4, extracted at r = 50 A1. Note the 
quasinormal ringdown exhibited by the linear decay in the 
amplitude and the constant frequency on this log-linear plot. 



if'^, a and /?' reads: 



(6) 



(7) 



where we have introduced the auxiliary variable ^E" in or- 
der to reduce the system to a first-order-in-time form. 
Notice that we label the time coordinate in the above 
two equations as t, in order to stress the distinction be- 
tween this and the time coordinate t associated with the 
BBH evolution. 

For each of the scalar evolutions, we provide initial 
data for and ^ as follows: 



4 - 
r e 



^' = - 



2(r-ro) 



(8) 
(9) 



where rg = 1.6A^ and a = 1A4. An outline of our initial 
setup is provided in figure [H 

Notice that, as observed in ^8|, the fact that 
Re[Yfm] = [Yim + Yi,-m)/'2 implies the concurrent pres- 
ence of the {i, to) and (£, — m) modes in the initial data. 
In order to produce ringing waveforms with a single dom- 
inant frequency, we chose only to = modes. Follow- 
ing Dorband et al. [3] , we extract the multipoles of the 
scalar field at a coordinate radius r = 5Ai in all cases 
(the effect of the extraction radius on the quasinormal 
frequency is illustrated in section BlI Cp . 



E. Extracting black hole properties from QNM 
ringing 

In the discussion above, we have mentioned how scalar 
perturbation theory on a black hole background predicts 
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FIG. 2: Initial isosurfaces for the scalar field $ as in equation 
(IHl), in the 1 = 2 case (broken lines). The solid outline rep- 
resents the apparent horizon's intersection with the xy-plaxie 
(top) and its intersection with the plane containing the 2-axis 
and the line connecting the two black hole centers (bottom), 
at t = 160 A1. Given the initial data for ^ in equation ([9}, 
the two lobes along the z-axis start moving towards the origin. 
The detector sphere used for extraction is located at r = 5. 



that, as soon as the initial impulse driving the perturba- 
tion propagates away, the scalar field will relax back to 
its equilibrium configuration through a period of damped 
(quasinormal) oscillations; numerical simulations (see, 
e.g., p^ ) confirm this picture. In the quasinormal regime 
the scalar field can be quite accurately described in terms 
of a set of quasinormal modes identified by the three 
numbers ^, m (labeling the angular eigenstates) and n 
(labeling the radial eigenstates). In time, each multipole 
evolves perturbatively, independently of the others and 
is characterized by a specific complex frequency r^^mn, 
whose dependence on the black hole mass and spin is 
known. In principle, then, it should be possible to extract 



information about the black hole parameters M and j by 
measuring ^imn for a number of different £mn-modes. 

The peculiar nature of quasinormal modes complicates 
this procedure: since quasinormal modes do not consti- 
tute a complete set of eigenstates, they cannot serve as 
a basis for the solutions of equation and there is no 
guarantee that, at any given instant of time, the scalar 
field can be described as a superposition of quasinormal 
modes alone. However, the practical evidence is that the 
system will follow a time evolution path which, for some 
transitory period, will live almost entirely on the sub- 
space spanned by the quasinormal modes [2^, [2^ . 

As in Q, we compute the time series ^^Qir): 



$ffi(T) 



Mp,d,(f>,T)d^n (10) 



(where Yeo are the standard (spin zero) spherical harmon- 
ics with m = 0, and the angular modes $fo £^re extracted 
on a sphere Sp of fixed coordinate radius p) , and fit it to 
a superposition of quasinormal modes: 



<i>f^{T;A,M,j,cl),To) 



Re 



-iOj.o„(Mj")(T-rn)+0{on 



(11) 



where A£q„ and 4>eon are real constants, and the com- 
plex quasinormal frequencies ^eon{M,j) = LOgon{M,j) + 
iaion{M,j) depend on M and j in a somewhat compli- 
cated way, which can however be represented quite con- 
veniently through the use of tabulated data |29|, com- 
plemented with third order interpolation. Note that we 
only fit to the m — fundamental mode (n = 0), how- 
ever, we study the error associated with not including 
higher overtones. 

In order to extract the black hole parameters we min- 
imize the residual of the fit 



q{AM,j,<j),To,Tf) = 

O dr I $,0 (r) - ^- ) I 



(12) 



with respect to A, M, j and (j) (with the M and j depen- 
dence given by the table interpolation method mentioned 
above). 

As mentioned above, there is no a priori prescription 
for the initial and final times tq and Tf delimiting the fit- 
ting window in which the minimum search is performed. 
In [18] and [9], tq is included as an additional fitting pa- 
rameter, along with some heuristic prescription for Tf. 
In [T(J| , the value of tq was obtained by cross-correlating 
different angular modes, with two different cutoff crite- 
ria for Tf. In this work, we start by choosing the ini- 
tial and final times tq and Tf so that the fitting window 
excludes the initial prompt response phase and the fi- 
nal wave portion contaminated by numerical and outer 
boundary error, or tq ^ 25 Ai and r/ ~ 55A4. We then 
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perform a two-dimensional analysis of the variations in 
the best-fit parameters induced by a change in tq and Tf . 
This analysis provides an estimate of the frequency ex- 
traction error and ensures that small variations on either 
of these parameters have no appreciable influence on the 
fit results. 

As noted in [42l|, the extraction of the complex expo- 
nentials from a quasinormal ringing waveform is a del- 
icate procedure. In order to minimize equation (jl2p . 
we have tested the non-linear least-squares Levenberg- 
Marquardt algorithm, which proved to be inadequate 
for our purposes. We then chose a fitting algorithm 
based on an optimization method known as the Covariant 
Matrix Adaptation Evolutionary Strategy (CMA-ES), 
whose details can be found in [43] and which is avail- 
able as a convenient open-source Matlab routine [3] ■ As 
a consistency check, we have tested the fitting procedure 
employing CMA-ES, the Matlab Levenberg-Marquardt 
nlinf it routine, and obtained comparable results when- 
ever the latter routine converged to a minimum. We 
have also tested CMA-ES on synthetic waveforms such 
as those discussed in [Z^l , obtaining similar results to the 
Kumaresan- Tufts and matrix pencil methods presented 
there. 

Notice that equation pT|) involves an infinite set of 
overtones, i.e. different n-modes. In practice, however, 
the excitation of the vast majority of modes will be highly 
suppressed, and we will only have to consider the over- 
tones with contributions above the level of numerical er- 
ror. Recent evidence [9J shows that, even though higher 
overtones are usually heavily suppressed, their inclusion 
in the fitting routine noticeably improves the quality of 
the fit. Since our primary interest is in extraction of the 
black hole parameters, however, we have decided to focus 
on the frequency of the fundamental mode alone, choos- 
ing a To sufficiently large to effectively guarantee that 
the contribution of the overtones (which damp quickly) 
is below our numerical error. We discuss the error contri- 
bution from using the fundamental mode below in section 

Unci 

Other important error sources include the finite radial 
and angular resolution, the presence of numerical error 
from finite differencing, and possibly additional effects 
due to the dynamical slicing of the spacetime which ordi- 
narily occurs during a 3-1-1 evolution. In order to address 
these potential issues, in section fill CI we try to quantify 
their effect. The mass and spin estimates from indepen- 
dent sources also provide evidence that our cumulative 
error budget is correctly estimated. 



III. RESULTS 

As a first test of the validity of our setup we evolved 
a single Schwarzschild black hole in isotropic coordi- 
nates using the standard moving puncture recipe that 
leads to coordinate dynamics. We then took snapshots 
of the gravitational field variables and solved the scalar 



wave equation numerically on these hypersurfaces and ex- 
tracted the fundamental quasinormal frequency found in 
the scattered scalar field for a variety of choices of tq and 
Tf in Equation (jlip . The values of the real and imaginary 
part of the quasinormal frequency coincided with the re- 
sult from perturbation theory up to Aui < 0.002A^~^ and 
Aa < 0.002A1"^ (which corresponds to AM < O.OIM 
and Aj < 0.05 for the £ = 2 mode). 

In the binary coalescence runs, for all t > 160A1 a 
common apparent horizon is present, and becomes ap- 
proximately axisymmetric at i ~ 165M (see '451 fo'" ^ 
description of the Killing vector field finding algorithm 
we deployed). Starting at this time, and until after the 
coalescence is complete (at about t = 2607W , as indicated 
by the horizon parameters M and j becoming constant 
in time), we freeze the gravitational variables at regular 
intervals and perform an evolution of the scalar field as 
described by equation ([T]) on the corresponding geometry, 
using the reference resolution of A^/51.2. 



A. Waveforms 

Figure [3] shows the i = {2,4,6} modes of the scalar 
field extracted on a coordinate sphere of radius p = 5A4 
in separate panels. Each mode was evolved on three dif- 
ferent hypersurfaces labeled by t = {160, 200, 260}A1 (as 
stated above, at t — 160A1 the first common apparent 
horizon is found, and by t — 26QM the system has set- 
tled down and we expect to observe the typical evolution 
of a scalar field scattered from a Kerr black hole) . As can 
be clearly seen from this figure each waveform undergoes 
a damped oscillatory phase. We observe the emergence 
of quasinormal ringing on each spatial slice containing 
a common apparent horizon, regardless of whether the 
merged system has already settled down {t — 2607W) or 
is still evolving very dynamically (t — 160A1). Initially 
the scalar field ringing is virtually identical on each of 
the three hypersurfaces shown, but significant dephasing 
is obtained towards the end of the scalar field simulation 
due to the different values of the fundamental frequency 
in the three cases. This is not an artifact of resolution, 
but rather shows the difference between the hypersur- 
faces. 

For each time t where we extract a hypersurface snap- 
shot and evolve the scalar field, we fit a quasinormal 
mode to the scattered scalar field. Figure |4] shows the fit 
residual g as a function of t for the £ = {2,4,6} modes. 
This gives us a rough indicator of how well the scatter- 
ing can be described by quasinormal mode ringing. The 
fit residual q decreases with increasing t, indicating that 
the evolution on earlier hypersurfaces leads to an inferior 
fitting quality. This loss in fit quality at early times is 
reasonable because the scalar field is expected to depart 
from the familiar quasinormal ringing behavior since ear- 
lier hypersurfaces are representing a black hole in more 
and more dynamical phases. Also note that the higher 
modes are more difficult to resolve and hence a larger 



6 



10^ 



O 



o 

to 



10^ 




10^3 
10^ 

10-1 
10-3 



t/TW = 160 
t/M = 200 
t/M = 260 



VA^ = 160 
t/M = 200 
t/A^ = 260 




t/M = 160 
t/M = 200 
t/AI = 260 




FIG. 3: The £ = {2, 4, 6}-modes of the scalar field are shown in the three panels. Each mode was evolved on hypersurfaces 
extracted at t = {160, 200, 260}A1 from a binary black hole inspiral. The modes were obtained on a sphere of constant 
coordinate radius p = 
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FIG. 4: The minimum of the fitting residual q (see Equa- 
tion (|12p ) as a function of t for the i = {2, 4, 6} mode of the 
scalar field. At early times t the black hole has not yet set- 
tled down and hence the fit residual q is larger as the scalar 
field evolved on the black hole background does not exhibit 
similarly clean ringdown as it does for later times. 



residual is observed. 



B. Black hole parameters extracted from the scalar 
field probe 

Figure [5] shows the black hole mass M and spin j as 
obtained from the scalar field evolved on different hyper- 
surfaces labeled by t, for the i = {2,4,6} modes. For 
reference, we also show the values obtained directly from 
the dynamical horizon finder. Both M and j show quali- 
tatively similar behavior, with an initial transition regime 
immediately after the first common apparent horizon for- 
mation (at t = 160A^). About 20A^ after the common 
apparent horizon is found we can reliably extract the spin 
and mass of the final black hole from the scattered scalar 
field. 

We want to emphasize that, while the scalar field 
modes on the different hypersurfaces appear almost iden- 
tical, the extracted mass and spins show much larger dif- 
ferences due to the strong dependence of AI,j on the 
extracted quasinormal frequency flim (notice that the 
error bars in figure [5] represent not only the uncertainty 
associated with the best-fit procedure, but also the con- 
tribution to the error due to finite numerical resolution, 
to angular mode extraction at a finite distance from the 
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black hole and to the choice to neglect any overtones 
or different angular modes in the frequency extraction; 
these contributions will be discussed in section ITlI Cp . 

In the transition regime between and 180A^ the 

fit errors (in particular, those from varying the choice of 
To and Tf) become so large that an accurate extraction 
of the black hole parameters is not possible anymore. 
Note that the extraction of M and j at early times is 
so sensitive to tq and Tf that the associated uncertainty 
dominates the bars shown in figure [5] for t < 180A^. 

We conclude that our scattering procedure constitutes 
an accurate probe of the final black hole potential (de- 
fined in the perturbative sense as summarized in section 
ITT|) at times t ~ ISOAl, consistent with the results from 
the isolated horizon parameters. Between IGOAl and 
180A^, the black hole parameters M and j cannot be 
reliably extracted; however, the scalar wave is still show- 
ing quasinormal ringing. In this regime, due to the level 
of uncertainty associated with our numerical resolution 
and frequency extraction routine, we can only provide a 
broad constraint for the extracted parameters. 



C. Effect of finite resolution, extraction radius and 
number of fitting modes 

At this point we would like to discuss the different error 
sources for the final extraction of the black hole parame- 
ters M and j from the scalar field waveforms. We discuss 
how sensitive the extraction of these parameters is to the 
fit range [to,tj], to resolution, to wave extraction radius 
and to fitting the waveform to a superposition of more 
than one quasinormal mode. It is important to keep in 
mind that the black hole mass and spin parameters show 
a strong dependence on the scalar quasinormal frequen- 
cies 51 and hence a small fractional error in f2 translates 
into a much amplified fractional error in AI and, espe- 
cially, j. 

In order to provide a description of the fitting error, 
especially due to the choice of the fitting window [tq, r/], 
we select tq = 25Ai and Tf = 55Ai as reference values 
and vary these two parameters over a lOA^-wide interval 
centered around each. The induced variation in M and 
j is illustrated in figure [S] for a few representative cases: 
this estimate provides a quantitative range for the un- 
certainty in the fitting parameters. In the following, we 
limit ourselves to this pair of reference values and quote 
the error associated with this choice. 

For the three different resolutions listed above {h = 
M/U.8, M/51.2 and M/57.6), we find that both the 
real and the imaginary part of the extracted frequen- 
cies exhibit a monotonic trend as h decreases, with the 
difference between neighboring resolutions decreasing as 
h ^ and the difference between the highest and the 
lowest resolution remaining always under 0.0003A^~^. 
As for the extraction radius, the six different choices {p = 
10M,20M, . . . , 50M in addition to the initial p = 5M) 
lead to frequency shifts of the order of 0.005A4~^. The 



uncertainty is also smaller if the extraction is performed 
on detectors that are farther from outer and refinement 
boundaries, such as the ones at 5A4, lOM and 30A^ (the 
different spatial resolutions characteristic of different ex- 
traction radii are not dramatically relevant, indicating 
that the modes are well resolved in all the cases consid- 
ered). 

Finally, the actual eigenfunctions of the angular part of 
the wave operator in equation ([T]) on a Kerr spacetime are 
the (spin-zero) spheroidal harmonics Sim-, rather than the 
spherical harmonics Ygm used here; however, as discussed 
inM,y^m = Sem + {a^m^)j:i'^iBi',nSi>,n+Oia^m, 
so that the use of spherical harmonics only causes a mod- 
est amount of mode mixing in the extracted waveforms 
(additionally, Berti et al. 147] find that this expansion is 
surprisingly accurate out to values of a|r2| close to unity). 
To determine the bias caused by such mode mixing on the 
best-fit value of the fundamental frequency, we extend 
the fitting function (jlip to a superposition of two and 
three different angular modes. The change in the com- 
plex frequency of the fundamental mode due to the inclu- 
sion of a different number of modes is under 0.0006A^~^. 
Unlike the error due to the frequency fitting, these last 
three sources of error have virtually no fluctuations from 
hypersurface to hypersurface. For all times, the cumula- 
tive errors associated with the determination of the com- 
plex frequency (including the last three sources of error 
but not the uncertainty introduced by the fitting pro- 
cedure) then amount to ~ 0.005A1~^. This number is 
propagated to M and j via AM = M{n + AD.) - M{n) 
and Aj ~ j{n + Ail) — j{n) which yields: 

AM = 0.02M , Aj = 0.1 for £ ^ 2 
AM = 0.02M , Aj = 0.15 for £ 4 
AM ^ OMM , Aj = 0.2 for £^6 . 

Finally, these numbers are summed in quadrature to the 
errors in M and j due to the fitting procedure (shown in 
figure [6l), yielding the error bars in figure [5] 

IV. DISCUSSION AND CONCLUSIONS 

In the previous section, we have presented the result 
of an experiment involving the scattering of a massless 
scalar field off the curvature of a spacetime containing 
two coalescing black holes. This spacetime was generated 
numerically by a full 3D black hole simulation following 
the moving puncture paradigm. Each spatial hypersur- 
face obtained served as a fixed background for the prop- 
agation of the scalar field as described in section lll Bl and 
the behavior of the scalar field was analyzed in each case. 

We observed the emergence of quasinormal ringing on 
each spatial slice containing a common apparent horizon 
regardless of whether the merged system had already set- 
tled down to a single Kerr black hole. While we could not 
reliably extract the final black hole parameters immedi- 
ately after the formation of a common apparent horizon 
(t = 160A^), we could still clearly identify quasinormal 
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FIG. 5: The mass and spin parameters extracted from the fundamental mode frequency, for the three angular modes. The 
error bars for the horizon mass and spin is included in the curve width. At early times, different choices of ro and r/ lead to 
significant modifications in the behavior of M and j as a function of t. The data error bars in this regime should therefore be 
considered as an indication of the mass and spin range on each hypersurface. 



ringing even in this early phase where the merged black 
hole is strongly excited and the scattering potential de- 
viates from that of Kerr. This finding agrees with earlier 
results where scalar perturbations on Vaidya spacetimes 
were studied and quasinormal ringing could be identified 
even on time-dependent black hole backgrounds [l^ . Our 
results and interpretation are also consistent with the ev- 
idence from the Close Limit approximation [ll|, where 
a pair of black holes in a head-on collision was shown 
to have already entered the perturbative regime when a 
common apparent horizon surrounded them. Based on 
this observation, the scattering of a test scalar field from 
the merged black hole could provide information about 
the final black hole state before the spacetime has settled 
down to its final state. 

Soon after horizon formation (i.e. for times t ^ 180A1) 
we were able to obtain estimates for the mass and spin 
of the final black hole from the scalar field probe in good 
agreement with direct isolated horizon measures. This 
indicates that, at this time, the dynamics surrounding 
the black hole has settled down sufficiently that the scat- 
tering potential for the scalar field is essentially that of 
the final black hole, even though there is still non-trivial 



dynamics in the vicinity of the horizon. 

A more exhaustive analysis of the early-time scalar 
quasinormal ringing waveforms before t = 180A^ is nec- 
essary to obtain more conclusive results in the non- linear, 
merger phase, and is left for future work. 
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